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1 INTRODUCTION 



ABSTRACT 

The spheroid of the Sombrero galaxy, NGC 4594, is considered a prototype of classical, 
merger-built bulges. We use a Spitzer, IRAC 3.6/im image to perform a detailed struc- 
tural analysis of this galaxy. If one fits to this image only bulge and disc components, 
the bulge occupies a locus in the mass-size relation close to that of elliptical galaxies. 
When an outer stellar spheroid is added to improve the fit, the bulge Sersic index 
drops by a factor of ^ 2, and, if taken at face value, could mean that this bulge is 
actually a disc-like, pseudo-bulge, or a bar viewed end-on. The bulge effective radius 
and the bulge-to-total ratio also drop dramatically, putting the bulge in a position 
closer to that of bulges in the mass-size relation. We discuss implications from these 
findings, including the locus of the Sombrero bulge in the black hole mass vs. bulge 
mass relation. With this new bulge mass estimate, current dynamical estimates for 
the mass of the central black hole in Sombrero are more than 10 times larger than ex- 
pected, if only the bulge mass is considered. A better agreement is found if the sum of 
bulge and outer spheroid masses is considered. Furthermore, residual images show the 
presence of a stellar ring and a stellar, inner ring or disc, with unprecedented clarity. 
We also show that Sombrero is an outlier in scaling relations of disc galaxies involving 
the disc, the spheroid and the globular cluster system, but not so when its structural 
components are considered independently. In this context, the globular cluster system 
of Sombrero might not be representative of disc galaxies. Finally, we discuss the possi- 
bility that Sombrero formed as an elliptical galaxy but accreted a massive disc, which 
itself has secularly evolved, resulting in a complex and peculiar system. 

Key words: galaxies: bulges - galaxies: evolution - galaxies: formation - galaxies: 
haloes - galaxies: individual: NGC 4594 - galaxies: structure 



The Sombrero galaxy (NGC 4594, or M 104) is a lo- 
cal, massive disc galaxy. Classified as an unbarred Sa 
(de Vaucouleurs et al.lll99"ll V it is at a distance of 9.1Mpc 
(Jensen et al. 2 0031 , from infrared surface brightness fluc- 
tuations), and has a stellar mass of - 2.3 X lO^^M© (e.g. 
iTempel k, Teniej|2006t V At this distance, 1" corresponds to 
about 44pc. It is seen al most perfectly edge- o n, with an in- 
clination angle of ~ 84° (Emse llem et al.lll99 6). It exhibits a 
remarkable spheroidal structure which extends to distances 
much larger than usually seen in other inclined disc galaxies. 

Figure [T] shows an archive Spitzer IRAC image of this 
galaxy at 3.6/im ffrom iKennicutt et al. I [20031 ). where there 
is very little dust absorption or emission, as well a s little 
contamination from hot, young stars (see e.g. She th et al.l 
I2OIOI ). It is thus a superb source to study the bulk struc- 
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tural properties of the galaxy. This deep image shows that 
the extended stellar halo seems rounder than the more com- 
pact central bulge, closer to the disc plane. In early stud- 
ies based on shallower or optical images the full extent of 
the halo is not evident, and thus Sombrero has been tradi- 
tionally considered a bulge+disc systen i only, wit h a very 
high bulge-to-tota l ratio {B/T = 0.86 in iKentll 19881. 



see also 



iLarsen et af]l200ll : IPrake et al.ll2003l : iGiiltekin et al. 2009*) 
This difference between the shape of the extended halo and 
the central bulge might suggest different formation processes 
for both components. Alternatively, if halo and bulge form 
as a single entity, then they could have followed different 
evolutionary paths. A further alternative i s that the va ria- 
tion in ellipticity is a projection effect fsee lRvdenlll99lh on 
a single spheroidal component. 

Since most previous works did not distinguish the halo 
from the bulge, the whole spheroid was well fitted by a Sersic 
law with high Sersic index, i.e. n ^ 4, although then the 
disc had to be fitted using complicated non-standard func- 
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2 STRUCTURAL ANALYSIS 



2.1 Fits to spheroid and disc 




Figure 1. Top: IRAC channel 1 image of the Sombrero galaxy 
(NGC 4594) at 3.6/im, used in this study. Bottom: isophotal con- 
tours from the same image. Foreground stars are removed. The 
stellar halo stands out clearly and appears rounder than the cen- 
tral bulge. 



tions fsee iKentI [iQiil l . The spheroid is obviously rounder 
and vertically more extended than the disc, and it also 
has no conspicuous star formation. For these reasons, the 
spheroid in NGC 4594 has been considered a prototy pe of 
classical bulges (e.g. in iKormendv KennicuttI l2004l l , i.e. 
those formed through violent processes, as opposed to disc- 
like bulges, which are thought to form from disc material 
brought to the centre via non-axyssimetrical components 
(suc h as bars), which f orm as a result of disc instabilities 
(see Athanassoulal 1 20051) . In fact, it has been advocated by 
e.g. iFish er Drorv (2008) that the bulge Sersic index can 
be used to separate classical bulges fr o m di sc- like bulges, 
with a threshold at n 2. In iGadottil (|2009l ), it is argued 
that, although a very useful parameter, n has to be used 
with care when distinguishing the different bulge categories, 
since uncertainties in its measurement can be relatively high. 
With a different methodology to separate disc-like bulges, it 
is shown that the median value of n for disc-like bulges is 
1.5, with a standard deviation of 0.9. For classical bulges, 
one finds that the median n — 3.4, with standard deviation 
of 1.3. Furthermore, disc- like bulges have typical values of 
B/T below 0.1, whereas classical bulges have a distribution 
of B/T that peaks at ^ OA. 

In this study, we perform a detailed structural analysis 
of the Sombrero galaxy, using the Spitzer image shown in 
Fig. [U in Sect. [21 We show that if an outer stellar spheroid, 
or halo, is included, both the bulge n and B/T drop dra- 
matically. We discuss the implications of this finding in Sect. 

n 



We used budda to decompose NGC 4594 into its differ- 
ent structural components. Firstly, we fitted a model with 
bulge and disc components only (model bd). The code uses 
the entire galaxy image and each component is described as 
a series of concentric generalised ellipses. The bulge follows 
a Sersic profile and 5 parameters are left free to fit it: the 
Sersic index, effective radius, effective surface brightness, po- 
sition angle and elliptici ty. The disc is fitted as an edge-on 
disc, following Eq. (7) in Ivan der Kruit Searlel (jigSlI ). and 
3 parameters are fitted: scale- length, scale-height and cen- 
tral surface brightness. Although the disc is not perfectly 
edge-on, tests showed that an edge-on disc model produces 
marginally better results in this case. The reader is referred 
to Ide Souza et al.1 (|2004l ) and IGadottil ()2008l ) for details on 
the workings of the BUDDA code. 

Two points deserve further discussion here. Firstly, 
the Somb rero galax y has a weak AGN, more specifically a 
LINER (Bendo et al . 2006). These authors produced 2D fits 
to images of Sombrero at wavelengths ranging from 3.6/im 
(using the same image as in here) to 850/im, and found that 
this nucleus is brighter at 24/im. At 3.6/im, due to prob- 
lems in modeling the PSF, they excluded fro m the fit the 
inner 10", thus excluding the AGN component. I Jardel et al.l 
(|201lD recently performed a bulge/disc decomposition of 
Sombrero and avoided the AGN by excluding the inner 0.17" 
from their fit. A discussion on how the structural analysis 
presented here compares with these studies is given further 
below. Typically, only bright type 1 AGN s hould be mod - 
eled in fits such as the one presented here fsee lGadottill2008l ) , 
and the effect of not including an AGN component in these 
cases is usually to overestimate the bulge-to-total ratio, as 
well as the bulge Sersic index. In order to check whether an 
AGN component should be included in the fits to the Spitzer 
3.6/im image, we produced a test fit with such component. 
As expected, since the AGN is faint, adding an AGN com- 
ponent yields very similar results as when the component is 
absent, and thus its inclusion is not relevant. The most sub- 
stantial change is in the bulge Sersic index, which is ^ 10% 
lower when the AGN component is added. Such a difference 
is however within the la uncertainty, and the fit produced 
is slightly worse. In this context, it is important to note that 
the spatial resolution of the Spitzer image is relatively low, 
with a p ixel size of 0.75" and a PSF FWHM of ^2'\ In 
IGadottil ([2008), it is shown that at low spatial resolution 
even the light from bright type 1 AGN can be smeared out, 
and an AGN component should not be included in this case. 
A further reason to not include an AGN component in the 
fits here is that it can artificially reinforce some of the re- 
sults presented below, namely that the Sersic index of the 
Sombrero bulge can be lower than expected and that a fit 
with bulge and disc only is not a good fit. Thus, the fits 
presented here do not include an AGN component. 

Secondly, a good modeling of the PSF is critical for 
this study, since if it is not well modeled the resulting bulge 
parameters from the fits can be w rong. The PSF is mo deled 
as a circular Moffat function (see iTruiillo et al. lliooil). The 
FHWM and [3 Moffat parameters were estimated from fits 
to several point sources in the field of the Spitzer image. 
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Figure 2. Left: surface brightness radial profiles (in arbitrary 
units) obtained through ellipse fits to the IRAC image and to 
the BUDDA model images of a decomposition including bulge and 
disc only. Right: the same radial profiles, but obtained through a 
cut along the disc major axis. Lower panels show residual profiles 
after subtracting the full BUDDA model profile from the galaxy 
profile. The two circles mark the positions of the breaks in the 
galaxy radial surface brightness profile discussed in the text. The 
red solid lines have the same meaning as for the solid black lines, 
but when the fit is done masking out the outer ring only. 

using standard routines in IRAF0. We have also produced 
test fits using the PSF image built for the Spitzer Survey 
of Stellar Structure in Galaxies (S^G - ISheth et alJ '2010) 
by combining stars in many Spitzer 3.6/xm images, which 
returned similar results. 

For reasons that will be clear shortly, we also produced 
fits including an exponential halo (model bdh) and a Sersic 
halo with Sersic index constrained to be larger than 2 (model 
bdh2). The other free halo parameters are: effective radius, 
effective surface brightness and ellipticity (the position angle 
is the same as for the bulge). In addition, a fourth model 
was produced, with bulge and disc components only, but in 
which the bulge ellipticity varies with radius (model vebd). 
The results from the fits are shown in Table [1] and Figs. [2] 
tod 

The left panel in Fig. [2] shows surface brightness radial 
profiles of NGC 4594 and model bd, obtained through ellipse 
fits to the isophotes of both images, using the ellipse task 
in IRAF. If considered alone, it indicates that a good fit can 
be obtained with bulge and disc components only. 

However, in the case of edge-on galaxies, using ellipse 
fits to study the light distribution in galaxies can be mis- 
leading. In most galactocentric distances, the isophotes in an 
edge-on galaxy can result from more than one component. In 
the case of Sombrero, it is easy to see that, except from out- 
ermost isophotes, many isophotes are drawn from light com- 
ing from the disc, bulge and halo (see Fig. [1]). Thus, in this 
case, surface brightness radial profiles from ellipse fits are 
unsuitable to assess how light is distributed amongst differ- 
ent structural components, and hence to carefully judge the 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 3. Same as Fig. [2] but for a BUDDA decomposition includ- 
ing bulge, disc and an exponential halo. 
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Figure 4. Same as Figs. [2] and [3] but for a BUDDA decompo- 
sition including bulge, disc and a Sersic halo with Sersic index 
constrained to be larger than 2. 

quality of a 2D fit, since at each given radius such profiles re- 
spond in a complex fashion to how the different components 
influence the isophotes. Note that in the case of face-on or 
less inclined galaxies, although at any given isophote there is 
likely a contribution from more than one structural compo- 
nent only, such components are usually strongly dominant at 
specific radii intervals, which overcomes this difficulty, and 
render surface brightness radial profiles from ellipse fits to 
isophotes suitable. 

Thus, in order to more properly evaluate how light is 
distributed amongst the different structural components in 
the case of Sombrero, we made cuts along the major and 
minor axes of the disc to produce another set of surface 
brightness radial profiles (see right panel of Fig. [2] - the 
minor axis cuts are discussed in Sect. 12. 2|) . In contrast to 
profiles built via isophotes, in this second set of profiles, 
each radius correspond essentially to a unique galactocentric 
radius. Thus, in this case, each component dominates over 
the others at different radii intervals. Now one can see that 
model BD is a relatively poor fit to the outer parts of NGC 
4594. 

In the right panel of Fig. [2l the galaxy profile exhibits 
two important breaks (indicated by circles): the first, at 
r ~ 28", marks the limit inward to which the bulge domi- 
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Table 1. Results from the decompositions. Model BD has bulge and disc only. Model 
BDH includes an exponential halo, while model bdh2 includes a Sersic halo with Sersic 
index constrained to be larger than 2. Finally, model vebd has a bulge with varying 
ellipticity and a disc. Models sl to s3 are discussed in Sect. 12.2] and are fits to the galaxy 
spheroid only. Model Sl is a single Sersic spheroid, model S2 is a single Sersic spheroid 
with ellipticity fixed at 0.1, and model S3 is a single Sersic spheroid with varying ellipticity 
as in model vebd. 
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0.23 
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41 
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0.61 
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4.3±0.5 


69 
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0.32 




5.5 


Sl 






3.3±1.5 
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0.38 














S2 
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2.5±0.8 


81 

















Column (1) shows the model fitted. Columns (2) and (3) are the disc scale-lenght and 
scale-height, respectively. Columns (4) to (6) show the bulge Sersic index, effective 
radius and ellipticity, respectively. Columns (7) and (8) display the halo effective radius 
and elliptictity, respectively. Columns (9) to (11) are, respectively, bulge-to-total ratio, 
disc-to-total ratio and halo-to-total ratio. Finally, column (12) is the reduced of the 
fit. All scales are in arcseconds and the uncertainties in the bulge Sersic index are also 
shown. Check how the bulge ellipticity varies with radius in model vebd at the middle 
panel of Fig. (6] 



nates over the disc, while the second at r ~ 215", marks the 
limit outward to which the outer spheroid dominates over 
the disc. Note that the alternative interpretation in which 
the second break indicates where the bulge starts to domi- 
nate again over the disc is weakened by the fact that model 
BD does not fit the outer profile well, even with a bulge with 
Sersic index n = 3.9. This is due to the fact that the outer 
profile can be very well described by an exponential func- 
tion. Between these two breaks is where the disc dominates, 
and, in this region, one can see two bumps that are likely 
caused by two different stellar rings, which will be discussed 
further below. In addition, the middle panel in Fig. [6] shows 
that the ellipticity of the outer galaxy isophotes is substan- 
tially lower than those of the bulge component in model 
BD, which has an ellipticity e = 0.42. Furthermore, the cor- 
responding residual image in Fig. [71 shows clearly that the 
bulge component in model BD is also too luminous in the 
inner parts. Note that fitting an image where the outer ring 
is masked out produces very similar results (see the red lines 
in the right panel of Fig. [2]). This test was done by mask- 
ing only the outer ring component, and is different from the 
test performed in Sect. 12.21 in which the entire disc, ring 
and central sub-structures are masked out. In a ID, profile 
fit, the ring could produce a systematic effect, but because 
Sombrero is nearly edge-on, and the fits shown here are 2D 
fits, such systematic effect is absent, as the ring occupies a 
very small fraction of the pixels in the whole image being 
used by the code to produce the fit. It is also worth noting 
that BUDDA produces fits to the average, underlying distri- 
bution of li^ht from each com ponent (see e.g. Gadott Fet al.l 
l201Ql : [ScannaDieco et al. I I2OIQI ). Further, its minimization 
algorithm does search for the global x^ minimum through 
simulated annealing, avoiding to be trapped in a local x^ 
minimum. In the Appendix, we show that including two ring 
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Figure 5. Same as Figs. [2] to |4] but for a budda decomposition 
including a bulge, with ellipticity varying with radius^ and a disc, 
but with no halo. 



components in the modeling does not lead to significantly 
different results. 

We thus added an exponential halcEl to produce model 
BDH (Fig. [3]). The ellipse fits analysis suggests in this case 
an equally good fit, as compared to model BD. However, the 
cuts along the major axis, which, as argued above, is the 
proper tool to use in the case of Sombrero (together with 
cuts along the minor axis - see below), indicate that model 
BDH is an improved fit. The maximum in the residual profile, 
excluding the rings, drops by a factor of ~ 2. The reduced 



^ We use the word 'halo' to describe an additional component, 
whose presence is justified from a structural viewpoint. We cau- 
tion the reader that this does not necessarily imply the existence 
of a distinct constituent in terms of stellar population content 
(e.g. stellar ages and metallicities). 
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Figure 6. Radial profiles of position angle (top), ellipticity (cen- 
tre) and Fourier b4 coefficient obtained through ellipse fits to the 
IRAC image (points with error bars) and to the 3 full budda 
models, as indicated. There is practically no difference between 
models bdh and bdh2 concerning these profiles. The middle panel 
also shows the ellipticity profiles of the bulge component only and 
of the total model in model vebd. 



also drops by a factor of two. Note that this improvement 
of a factor 2 happens when moving from 10 free parameters 
to 13 free parameters, i.e. an increase of 30% in the number 
of free parameters. 

For the sake of exploring different alternatives, in model 
bdh2 we tried a Sersic halo, forcing a Sersic index larger 
than 2, as otherwise it would drop to about 1, as in model 
BDH. The best fit Sersic index for the halo in this case is 2.4. 
Figure [H shows that the fit in this case is not good in the 
outer parts, as the halo is too luminous. 

Figure [6] shows clearly how adding a halo to the model 
results in a much better fit to the geometrical properties of 
NGC 4594, as compared to model bd. The model images in 
Fig. [71 show that model bdh reproduces better the top panel 
in Fig. (TJ One reaches a similar conclusion analysing the 
residual images in Fig. [71 The residual image from model 
BD clearly shows that the bulge model is too luminous in 
the centre when the halo is not accounted for. As discussed 
above, this cannot be an effect arising from not including 
an AGN component, as this component, if at all necessary, 
is too faint. It is also unlikely to be a problem of PSF fit- 
ting, since we tried two different PSF models, which yielded 
similar results, and the area where the residual caused by 
the bulge is excessive has a size of almost ~1', i.e several 
times larger than the PSF FWHM, which we measured as 
1.9". It is thus clear that model bdh is a better fit, even if 
not fitting completely the galaxy image, for not having ring 
components. The residual images in Fig. [71 show clearly a 
ring in the disc of NGC 4594. Furthermore, it also shows 
an elongated central sub-structure, which could be an in- 
ner ring or disc. This is diffe rent from t he re latively faint 
nuclear disc already found by iBurkheadI (| 19861 ) , with a size 



of ?^15", which can also be seen in our residual images, but 
does not produce a clear imprint in the radial surface bright- 
ness profile of the galaxy. Both such structures induce the 
bumps in the surface brightness radial profile of NGC 4594, 
seen for instance in the right panel of Fig. [3l in the region 
where the disc dominates the profile. It should be noted that 
these residual images have a very narrow display stretch, in 
order to emphasize discrepancies between the model fitted 
and the galaxy, as well as existing sub-structures. Further, it 
is worth noting that, because the fits shown here are 2D fits, 
where the whole image is taken into account, as opposed to 
ID profile fits, the residual images are also an appropriate 
tool to verify the goodness of the fits. 

The variation of ellipticity with radius in Sombrero, seen 
in the middle panel of Fig. [6] suggests another, simpler possi- 
bility. A model with bulge and disc only, but where the bulge 
ellipticity varies with radius, could perhaps produce a fit to 
the galaxy at least as good as model bdh. RvdenI (|l99lh 
studied the geometry of ellipsoids with ellipticity varying 
with radius in different projections. Her figures 3 and 4 show 
a radial behaviour of ellipticity that resembles what we ob- 
serve for Sombrero. We thus tailored budda to allow for a 
bulge with varying ellipticity e, in order to produce model 
VEBD. In this model, thus, as in model bd, the word 'bulge' 
refers to the whole spheroidal structure seen in the galaxy, 
from the central bulge to the outer halo. However, to not in- 
troduce an excess of degenerated solutions, the way e varies 
with radius is fixed. The ellipticity variation shown in Fig. 
[6l gives us information on how this should be fixed, but one 
has to keep in mind that the variation measured through 
ellipse fits corresponds to the whole galaxy, not the bulge 
component alone. In order to evaluate the intrinsic elliptic- 
ity variation of the bulge component alone, one can look at 
the outermost radii, where the contribution of the disc com- 
ponent drops substantially. The ellipticity profile in Fig. [6l 
shows that outside a radius of about 250 arcsec, where the 
outer spheroid dominates over the disc, e drops more or less 
linearly from ^ 0.33 to ~ 0.1. Inside this radius, the ellip- 
ticity of the isophotes reflects that of both disc and bulge. 
There is in fact almost a plateau in e around this radius of 
250 arcsec. A natural choice is thus to keep e = 0.33 for the 
bulge component inside this radius. This choice is also sup- 
ported by the fact the ellipticity of the innermost isophotes 
in the galaxy is 0.41, and therefore relatively close to 0.33. 
The green line in the middle panel of Fig. [6l shows how e 
varies with radius for the bulge component only in model 
VEBD. Note that this choice also leads to a better fit using 
model VEBD. If one leaves the bulge ellipticity rising linearly 
until the centre, the resulting model is too eccentric. What 
happens is that the disc, because it is almost edge-on, has an 
ellipticity of about O.Qj. The spheroid being rounder helps 
dilute this eccentricity, less so where the disc dominates. If 
the spheroid is also eccentric this dilution is reduced. In fact, 
one can see that even with a bulge model with a maximum 



^ Note that our disc model is viewed perfectly edge-on, and in 
this case the concept of ellipticity is not applicable as it is for face- 
on or inclined discs. The disc ellipticity is simply not a parameter 
in our models. The value of 0.9 is a refer ence value, and i t come s 
from the measure of disc inclination by lEmsellem et al.l (|l996l V 
The disc scale length and scale height determine how thin the 
disc is, and one sees that their ratio varies amongst our models. 
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Figure 7. Images of the full budda models (top) and the corresponding residual images (bottom) after subtraction of the models from 
the IRAC image. From left to right: bulge+disc model, bulge+disc and exponential halo model, bulge+disc and Sersic n > 2 halo model. 
The top images are in the same scale and display levels as in Fig. [1] The residual images are all in the same scale and display levels. The 
one at left shows clearly that the bulge model is too luminous when the stellar halo component is not included. The two panels in the 
last row are the corresponding images for model vebd, i.e. bulge with varying ellipticity+disc. 



ellipticity of 0.33, the dilution is too weak to reproduced the 
galaxy ellipticity profile (blue solid line in the middle panel 
of Fig. El). 

Figure [5] shows the results of the fit from model vebd. 
It shows that, although the fit looks better than the one 
from model BD, model bdh, with its n ~ 2 bulge and an 
exponential halo, still provides a better fit than model vebd, 
with a bulge with varying e and n ~ 4, and no halo. In fact, 
the reduced from model vebd is similar to that of model 
BD, i.e. a factor about 2 worse than model bdh. Further, 
although the ellipticity of the full model vebd reproduces 
the more or less linear fall with radius, it does so with an 
offset from the values of ellipticity measured in the galaxy 
(Fig. [6]). It might be surprising that model vebd does not do 
much better, given its degrees of freedom. The crucial factor 
here appears to be the shape of the profile. Given that the 
outer part of the galaxy profile is quite strictly exponential, 
the fitting of a single spheroid leading to a Sersic index of 
about 4 results in a bad fit. Further below, we will discuss 
the possibility that the n ~ 2 bulge in model BDH is actually 
a bar seen end-on. This bar might be the one component 
which is leading model bp to a relatively bad fit. Figure 1 in 
[Bureau Athanassouial (|2005l ) shows unambiguously that 
an end-on bar remarkably resembles a classical bulge from 
a purely morphological viewpoint. 

The most striking result from this analysis is thus that 
the bulge Sersic index of the Sombrero galaxy drops from 
3.9 to 1.9, when one does a fitting of the IRAC image in- 
cluding the stellar halo in the model. Concomitantly, B/T 



drops from 0.77 to 0.1^, while the bulge effective radius 
Te drops from 3.1kpc to 0.44kpc. These drops in B/T, n 
and Te make the bulge in Sombrero look much more typi- 
cal, since the average values for Sa galaxies are B/T — 0.31 
(with a Icr scatter lower limit o f 0.15), n = 2.6 and Ve — 
0.7kpc ("Graham k WorlevI l2008l . from i^^-band data - see 
also Gadotti 2009 ). In t his co ntext, it is worthy to point out 
that iTempel k Teniesl (|2006h found a value of B/T 0.15 
through sophisticated dynamical modeling of the mass dis- 
tribution in NGC 4594, including a halo component, in very 
good agreement with our model bdh. Note that the value 
obtained for the bulge effective radius when the halo is not 
included in the fit is almost a factor 5 larger than the typical 
value for other Sa galaxies. Another way of looking at this 
is put in Fig. [H It shows that the spheroid in NGC 4594, if 
treated as a single component, resembles, from a structural 
point of view, more an elliptical galaxy than the bulge of a 
disc galaxy . In this figure, we show the mass-size relations 
obtained in'Cadotti' ("2009^ for nearly 1000 elliptical galaxies, 
classical bulges and pseudo-bulges, drawn from a complete 
parent SDSS sample, with a typical redshift of 0.05. The 
spheroid in Sombrero as given by model BD is larger and 
more massive than any of the bulges found in that sample. 
In fact, it is as massive as the most massive elliptical galax- 



^ I Graham &: WorlevI (|2008l ) discussed how early studies forcing 
the bulge Sersic index at 4 resulted in overestimated bulge-to- 
total ratios. 
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Figure 8. Mass-size relations for el liptical g alaxie s, classical 
bulges and pseudo-bulges from data in lGadotTil (|2009l ). The posi- 
tion of the Sombrero bulge in both models BD and bdh is marked. 
The spheroid in Sombrero as given by model BD is larger and 
more massive than any of the bulges. In fact, it is as massive as 
the most massive elliptical galaxies, and falls near the mass-size 
relation of ellipticals, rather than that of (classical) bulges. On 
the other hand, the bulge obtained in model bdh, although some- 
what small for its mass, falls near the relation drawn for classical 
bulges, although the locus it occupies in this relation is also not 
unusual for a disc-like, pseudo-bulge. 



ies, and falls near the mass-size relation of ellipticals, rather 
than that of ( classical) bulges. In th is context, the scenario 
entertained bv lEmsellem et al] |l996l V in which the spheroid 
forms after the dissolution of a bar, is highly improbable, 
given that the typical bar is about two orders of magnitude 
less m assive than the spheroid in Sombrero (see e.g. GadottJ 
l201lh . On the other hand, the bulge obtained in model bdh, 
although somewhat small for its mass, falls near the relation 
drawn for classical bulges, although the locus it occupies in 
this relation is also not unusual for a disc-like, pseudo-bulge. 

In summary, our structural analysis shows that the 
most appropriate model for Sombrero is model bdh, with a 
spheroid comprising an inner bulge and an outer halo. Nev- 
ertheless, we cannot rule out the possibility that the spheroid 
is in fact an elliptical galaxy (as suggested in Fig. [8|), that, 
unlike most elliptical galaxies, displays intricate structural 
properties, due to the presence of a massive embedded disc, 
and the possibility that this disc have formed a bar and thus 
experienced effects from secular evolution. This opens up the 
possibility that the current view on the Sombrero spheroid 
as a single entity, a large classical merger-built bulge, is mis- 
taken. These possibilities will be discussed in more detail in 
the next section. 



2.2 Fits to spheroid only 

Another important way to address the nature of the spheroid 
in Sombrero is by looking at minor axis profiles, where the 
influence of the inner sub-structures, disc and rings is mini- 
mized. In Fig. [9] we plot the minor axis profiles of the galaxy 
and models BD, bdh and vebd. One sees that model bdh 
provides a better fit, although it does not fit well the outer 
parts. This is very likely because the halo has a fixed el- 
lipticity of 0.28, somewhat larger than the ellipticity of the 
outermost parts of the galaxy, which is around 0.1. Model 
VEBD provides a better fit than model BD, attesting the im- 
portance of having an ellipticity varying with radius. At the 
inner parts, models BD and BDH provide equally good fits, 
whereas model vebd is slightly off. 

In this context, we produced further fits by excluding 
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Figure 9. Left: minor axis surface brightness profiles of NGC 
4594 and 3 of the models fitted, as described. Right: Same profiles 
with with the radial axis in logarithmic scale, to emphasize the 
inner parts. 



minor axis 




250 
r (arcsec) 



10 

r (arcsec) 



Figure 10. Same as Fig. [9l but for the models with only a 
spheroid. Model Sl is a single Sersic spheroid. Model S2 is a single 
Sersic spheroid with ellipticity fixed at 0.1. Model s3 is a single 
Sersic spheroid with varying ellipticity as in model vebd. The 
vertical line marks one disc scale-height, within which the galaxy 
image was masked out. 



all sub-structures apart from the spheroid. This was done 
by masking a rectangular area centred at the galaxy center, 
and with dimensions 36" x 430" , along the disc major axis. 
The minor axis dimension is twice the disc scale height (as 
in model bd), and the major axis dimension was chosen to 
allow only for the part dominated by the spheroid, as seen 
in the galaxy profile (see e.g. Fig. [2]). 

The new fits include a Sersic spheroid only, with the 
Sersic index left free. Model sl has free ellipticity, model 
s2 has a fixed ellipticity of 0.1, i.e. reproducing the outer- 
most parts of the galaxy ellipticity radial profile (see Fig. 
[6|), and finally, model s3 has varying ellipticity, as in model 
VEBD. The resulting parameters are shown in Table [1] Mi- 
nor axis profiles from these fits are shown in Fig. 1101 These 
new models generally provide better fits to the outer parts 
of the galaxy than the models in the previous sub-section, 
but fail in the inner parts. This evidently reflects the design 
of this exercise and the influence the disc component has in 
the global fits of the previous sub-section. 

Model sl is the one that fails the most in the outermost 
parts. Model s2 is overall quite good but fails at interme- 
diate radii, i.e. from about 50" to 100". Finally, model s3 
can account for the change in ellipticity in the outermost 
parts, but generally does not result in a very good fit. This 
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exercise illustrates the fact that even excluding from the fit 
all sub- structures apart from the galaxy spheroid, and even 
allowing for varying ellipticity, Sombrero's spheroid is not 
well described by a single Sersic component. This justifies 
adding an extra component apart from the main bulge and 
disc components. Furthermore, the left panels in Figs. [9] and 
[10] clearly show that the outermost surface brightness profile 
of Sombrero follows very well an exponential profile, as the 
halo in model bdh, and thus cannot be well reproduced by 
a Sersic function with high Sersic index. In fact, as shown 
in Table [U model sl has n — 3.3, while models s2 and s3 
have n — 2.5. This suggests that the change of ellipticity is 
related to a flatter surface brightness profile, since models s2 
and s3 not only have lower Sersic indices, but also reproduce 
better the outer lower ellipticities seen in the galaxy. This 
is a further indication of the presence of an extra spheroidal 
component, with flatter mass distribution and lower ellip- 
ticity. 



3 THE INTRIGUING NATURE OF 
SOMBRERO'S SPHEROID 

3.1 Evidence from structural analysis 

The work done in the previous section is a compelling il- 
lustration of the fact that structural analysis of galaxies 
through image decomposition is a powerful tool, yet rather 
a complex affair. It demonstrates that results can signifi- 
cantly depend on details that show up only through better 
data and/or careful scrutiny. In fact, ID profile fits alone 
can also be misleading, as shown in Fig. 7 of FGrahaml (|2QQlh . 
Our structural analysis indicates that the spheroid in Som- 
brero comprises two main components: an inner bulge and 
an oute r halo. 

Ben do et al.l (|2Q06h presented a similar study as ours us- 
ing the same IRAC image. They did not include a halo com- 
ponent in their fits. However, their focus was on the dust dis- 
tribution in Sombrero, and so they concentrated their efforts 
on the 24/im image, where, as they stated, the "bulge almost 
disappears completely". Although they showed a residual 
image from the fit to the 24/im image, they did not show 
intensity radial profiles, making it somewhat more difficult 
to assess how good was their fit. In particular, they did not 
show a residual image from their model to the 3.6 /xm image, 
which basically traces only stars, and where the spheroid is 
very conspicuous. We also note that they fixed the Sersic 
index of their bulge model to n = 4. Their B/T for the 
3.6jtxm im age is 0.77, i.e. the same as we find in model bd. 
Ijardel et al.l (j201l|) also presented a similar study, putting 
together HST and ground-based images. They used a non- 
parametric form to fit the bulge profile, plus 3 discs, and 
found B/T — 0.73. They stated that their bulge profile 
"could in principle be a combination of a Sersic bulge plus 
exponential halo". 

If the low value of the bulge Sersic index in the Som- 
brero galaxy is correct - as the results in this study suggest 
since model bdh gives us the best fit - and if Sombrero has 
indeed a classical bulge, it is an example of how using the 
bulge Sersic index as a discriminator between classical and 
disc- like bulges can lead to wrong bulge classifications. This 
is only surprising because NGC 4594 is a prototype for clas- 
sical bulges. Alternatively, the lower value for n found here 



could mean that the Sombrero galaxy has actually a disc-like 
bulge, rather than a massive classical bulge, which would be 
truly surprising. 

Two of the previous arguments in favour of the Som- 
brero bulge being a classical bulge are now weakened. Its 
Sersic index is not as high as thought, and it does not ex- 
tend out from the plane of the disc as much as when the 
halo was considered part of the bulge (although even the 
bulge in model bdh is not as flat as the disc). On the other 
hand, the lack of substanti al star formation, as sho wn e.g. in 
G ALEX observations (see ' Gil de Paz et al.l l2007f ) , suggests 
that Sombrero has a classical bulge. 

Perhaps the most important criterion to identify a disc- 
like bulge comes from kinematics: in contrast to classical 
bulges, disc- like b ulges have significan t rot ational support. 
'Kormendvl (|l988h .lWagner et all (| 19891 ) and'Emsellem et al.1 
(1996) find strong rotational support in the inner regions 
of Sombrero. This would support the suggestion that the 
bulge in Sombrero is actually a disc-like bulge. In this case, 
the concept of the Sombrero bulge as a prototypical clas- 
sical bulge is a result of the presence of its extraordinary 
stellar halo, whose full extent can only now be properly ap- 
preciated. Without knowing the true geometrical properties 
and light profile of this halo, it was immediately conceived 
as a massive, classical bulge. If we follow this reasoning, we 
conclude that, basically, what has been called 'bulge' is ac- 
tually a massive stellar halo. This halo, in contrast to the 
bulge, does not contribute considerably to the light/mass 
distribution at the galaxy central region. 

The residual images in Fig. [71 clearly show a stellar ring, 
apart from another central stellar substructure, which could 
be a nuclear ring or disc. Suggesting evide nce of these com- 
ponents have been publish ed before (see IWainscoat et al.l 
Il990l : temselle m etHl Il996l l. but only with the IRAC im- 
age they stand out so clearly. An ou ter dust rin g sho ws up 
beautif ully in another Spitze r image (|Kennicutt et al.ll2003l , 
see also I Vlahakis et aLll2008l V Rin gs are usually, but n ot al- 
ways, associated to bars (see e.g. iGrouchv et aLllioiol V An 
end-on bar in Sombrero can explain the central kinematics 
( Emse llem et a l. 1996), and wo uld be a natural explanation 
for the presence of the ring(s). lEmsellem FerruitI ()2000h 
went further and suggested that a nuclear bar could also be 
present. One cannot currently rule out the possibility that 
the bulge component fitted above in model bdh is (at least 
partially) actually a bar seen end-on. This would be con- 
sistent with a photometric bulge with low Sersic index and 
yet no young stellar popula tion, as bars are usually p opu- 
lated with old stars (see e.g. iGadotti de Souz ^ l2006h . An 
end-on bar in Sombrero would also explain why the central 
isophotes are not as flat as the disc, as one would expect for 
disc-like bulges. 

Using the value of B/T — 0.13 found in model bdh, 
the total galaxy stellar mass and the black hole mass vs. 
bulge mass relation given in iHaring Rixl (|2004l ). one finds 
that the expected black hole mass is 4.1 X lO^M©. This 
is in strong cont r ast to dynamical measurements given by 
Kormendv et al. 1 ([l996l), which point out a black hole with 
mass about lO^Mc?) . More recen t results also point to large 
black hole masses. Ijardel et al ] (|201lh found a black hole 
mass of 6.6 X lO^M©. If one adds the halo luminosity to the 
bulge, the expected black hole mass rises up to 2.5 x lO^M©, 
which alleviates the discrepancy. If one uses the bulge mass 
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in model bd, the black hole mass derived is 3 x lO^M©. In 
this context, it is interesting to point out that the dynamical 
black hole mass measurement is also at odds with the black 
hole mass vs. velo city dispersion rela tion (see iBeifiori et al] 
20091. but see also [jardel et alJboilh . Using the results in 



Giiltekin et al. 120091), the expected black h ole mass from 
the la tter relation is 2.9 x lO^M©, although I Graham et al.l 
(|201lh argue for a value of 5 x lO^M©. It thus seems that 
either more sophisticated dynamical modeling is necessary 
to properly understand the central mass concentration in 
Sombrero, or it truly has an exceptionally massive black 
hole. Alternatively, the halo mass has also to be accounted 
for. 

These results urge one to question thus how the halo 
in Sombrero did form, and why such structure is seen in 
some galaxies but not in all. For instance, NGC 4565 is also 
an edge-on galaxy, which was observed at t he same wave- 
length also using Spitzer (jSheth et al]l2010l ). more deeply 
than the Sombrero image used here, and yet it does not 
show such an extended halo. Conv ersely, NGC 5866, also ob- 
served bv'Kennicutt et"al] (|2003l l does exhibit an extended 
halo. Early merger ev ents are a likely formation mechanism 
for the Sombrero halo (jEmsellem et al.lll996l ). Therefore, the 
properties of such haloes and their host galaxies, and how 
often they are present in galaxies, can provide constraints 
to physically characterise merger events, and elucidate how 
often and in which circumstances they occur. This is evi- 
dently an important issue for theories of structure growth 
in the universe. In this c ontext, it is w orthy to point out that 
ISeigar et all (l2007l ) and lPierini et al.1 ^08) found that the 
haloes in BCGs (Brightest Cluster Galaxies) are better de- 
scribed with exponential luminosity profiles than with Sersic 
profiles with n 4, as in the case of Sombrero's halo. This 
also lends support to the arguments above that the main 
spheroidal in Sombrero is indeed a massive stellar halo. 

Finally, one is left wondering on the impact of such ex- 
tended haloes in our view of face-on galaxies. Since they 
might not be readily discernible in such cases, results from 
image decomposition can be substantially altered if they 
are not included in the model fitted, as seen above. Anti- 
truncated outer profiles can some times be a signature of 
the presence of a massive s tellar halo, as poi nted out by 
lErwin et all (|2005l . see also lErwin et al1l2008l ). Note that 
such a feature is clearly seen in the outer parts of Som- 
brero's radial surface brightness profiles shown in Sect. [21 
The rightmost circle in the right panel of Fig. [2] marks the 
position where the anti-truncation begins. 

3.2 Evidence from further scaling relations 

As mentioned in the Introduction, NGC 4594 is often re- 
garded as a prototypical early-type disc galaxy. With its 
massive spheroid and the almost edge-on dis c, it is usu - 
ally classified as an Sa spiral (de Vau couleurs et al.lll99lh , 
although some aut hors prefer to consider it an SO (e.g., 
iRhode ZeDll2004l . RZ04 hereafter). The results presented 
in the previous sections, however, raise the intriguing possi- 
bility that the spheroid may not be a classical bulge, but to a 
large extent a stellar halo or an elliptical galaxy (see Fig. 8). 
The latter possibility could in principle be considered just a 
matter of semantics if one takes the view that bulges are sim- 
ply ellipticals surrounded by a prominent disc fe.g. jRenzinil 



ll999l : lKormendv k Kennicuttil2004l \ but Fig. [8] shows that, 
at least at the high mass end, bulges and ellipticals follow 
different scaling relations. At fixed stellar mass, ellipticals 
tend to have larger effective radii than bulges, and it is note- 
worthy that Sombrero's spheroid as a whole lies much closer 
to the sequence defined by the former. While this does not 
necessarily preclude that they might share similar formation 
mechanisms, it clearly indicates that at least their formation 
histories must differ. 

The resemblance of Sombrero's sphero id with ellipti- 
cal galaxies has of course been noted earlier. iHes Peletieil 
(1993) pointed out that the mild colour and (absorption) 
line strength gradients, and the central kinematics, indi- 
cate such similarity. Indeed, the red colours of the spheroid 
inner regions gradually turn bluer with radius, matching 
those of the metal-rich glob ular cluster system (GCS) at 
about 5kpc from the centre (jSpitler et al. l2006l V This be- 
ha viour appears to hold out to large galactocentric radii, 
as iMould Spitlerl (|2010 ll found that the metallicity distri- 
bution function (MDF) of stars in a field at about 18kpc 
from the centre peaks at the same value ([Fe/H] ^ —0.5) 
as the metal-rich GCS - potentially indicating a link be- 
tween the formation of both components. This peak is sig- 
nificantly more metal-rich than t he halo of the Milky Way 
analog NGC891 ([Fe/H] ^ -1: iReikuba et al.ll2009l ). but 
comparable to what is found in other elliptical galaxies (e.g. 
NGC 5128; Harris & Harris 2002). 

The GCS in Sombrero perhaps provides the most com- 
pelling evidence of its spheroid's peculiar nature. It has long 
been known that NGC 4594 possesses the largest GC popu- 
lation ever found in a disc galaxy. RZ04 uncovered a GCS 
exceeding Ngc = 1900 globular clusters that extends out to 
^ 50 kpc, with ^ 40% of them being metal-rich. Figure [TTI 
shows the GC mass specific frequency0 for a compilation 
of > 10^° Mq elliptical, lenticular and spiral gal axies 



in di fferent environments - from field /group s (jSpitler et al.l 
I2OO8I I to the Virgo cluster (I Peng et al.ll2008l ). With a factor 
~4 more GCs than M31, and inhabiting a low-density envi- 
ronment, it is clear that both the mass and Tn of Sombrero 
are more characteristic of elliptical galaxies than of lenticu- 
lar s or spirals. This plot also indicates that a major merger 
event between two M31-like spirals cannot reproduce the GC 
abundance of NGC 4594 - not even that of the old, metal- 
poor population, which differs by a factor of about 3.5 be- 
tween M31 and Sombrero. S pitler et al. (2006 ) pointed out 
that Sombrero has the largest bulge specific frequency of red 
GCs among all spirals and thus its value is more characteris- 
tic of more massive galaxies. Moreover, they discovered that 
metal-poor GCs follow a colour-magnitude relation, making 
NGC 4594 the first di sc galaxy where su ch trend was ob- 
served. More recently, iHarris et al.l (|2010l ) showed that this 
mass- metallicity relation scales as Z oc L*^"^, pretty much 
equivalent to what has been found in most massive ellipti- 
cals. 

In summary, the spheroid in Sombrero has properties 
that differ significantly from those of classical bulges, in- 
cluding its mass-size relation; the richness of its GCS (and 
especially the metal-rich population) ; and the high metallic- 
ity of the outer halo when compared to other disc galaxies. 



Tn = A^n./rM^/10^Mr7.): [z^f fc AshmaiJ (ll993l V 
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Figure 11. Globular cluster mass specific frequency of ellipti- 
cal, lenticular and spiral galaxies in different environments (see 
text for details). NGC4594 has the largest GCS of all known 
disc galaxies, and occupies a locus that is more characteristic of 
intermediate-mass ellipticals. The black solid circle indicates the 
corresponding value if the whole GCS were associated just to 
the spheroid - i.e., not taking the disc mass into account. 



It is also pertinent to ask ourselves whether the 
embedded disc (see Fig. 1) displays peculiar properties. 
ICatinella et~all (|20ld ) show that the neutral gas mass frac- 
tion in the GASS sample strongly decreases as a function 
of stellar mass, with ^ 2 x 10^^ M© galaxies having a 
weighted median Mi///M^ = 0.016, with a 0.386 dex scatter 
(see Fig. [121 top panel). NGC4594, with a stellar mass com- 
parable to this value b ut a much low er neutral gas content 
(Mi// 3.1 X 10^ Mr:.: lBaiaia et al.| [T984FI) deviates ?^5cr 
from this relation. Considering instead the 5 x 10^° disc 
mass would bring this ratio to within ^2a oi th e observed 
relatio n in the GASS sample. As discussed by iHau et aP 
("2008'), this massive stellar disc implies that its peak SFR 
must have been a factor 5-10 times higher in the past than 
the current 0.42 M© yr"^ (|Skibba et al.ll201lh . and the halt 
of substantial disc growth is suggestive of an advanced evo- 
lutionary stage. This is supported by the analysis of stellar 
populations in the outer disc, that indicates old ages (>8 
Gyr) and close to solar metallicities - s imilar to those of 
the inner spheroid (jEmsellem et al ] |l996l ). Figure [I2] (mid- 
dle panel) shows the relation between specific SFR and stel- 
lar mass for differen t galaxy types in the KINGFISH sample 
(jSkibba et al.l 1201 il l. NGC4594 again stands out as a clear 
outlier within spiral galaxies, but has rather normal values 
for ellipticals. Perhaps not surprisingly, the stellar disc oc- 
cupies the same region as the other spirals when considered 
an independent entity, and the same is true for the dust-to- 
stellar mass content (Fig.[T2l bottom panel). This suggests 
that Sombrero's disc is not intrinsically peculiar but, in- 
stead, it is the out-of-proportion spheroid what makes the 
galaxy an outlier in these scaling relations. 

^ All relevant quantities refer to our adopted D = 9.1 Mpc dis- 
tance. 
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Figure 12. Scaling relations of NGC 4594 as compared with those 
of other > 10^ Mq elliptical, spiral and irregular galax- 
ies. Top: pentagons show the neutral gas to stellar mass con- 
tent for galaxies of all morphological type s with Hi detections 
in the GASS survey (Catinella et aP l201oh . Dashed lines show 
the weighted median and 3 a limits when non-detections are also 
included. Middle: Specific star formatio n rate of galaxies in the 
KINGFISH sample (Skibba e t 311120111 ), with different symbols 
corresponding to different morphological types. Bottom: same as 
for the middle plot, but for the dust to stellar mass ratio. In 
all three panels it is obvious that NGC 4594 is a clear outlier of 
the relations followed by spirals, but not so if the stellar disc is 
considered as an independent entity (solid symbols). 



3.3 A scenario for Sombrero 

The emerging picture of Sombrero supports a scenario where 
most of the formation activity took place at an early epoch. 
The bulk of stellar mass settled in a massive and extended 
spheroid at z > 2. This spheroid exhibits a mild metal- 
licity gradient, from roughly solar in the inner regions to 
a peak [Fe/H] ^ —0.5 in the outskirts - the same values 
as for the metal-rich GCS. This latter value can be repro- 
duced, in chemical evolution models that include gas in- 
fall, with an initial burs t of gas accretion lasting < 0.2 Gyr 
(iMouid Sp 'ityiioT^. Disc formation followed quickly, 
leaving ample time for secular evolution to take place. Dy- 
namical resonances due to the formation of a bar could 
have shaped the complex structures revealed in the resid- 
ual images previously shown, potentially contributing to 



Does Sombrero have a (classical) bulge? 11 



the s pheroid structure in the inner regions (jEmsellem et alJ 
Il996l ). 

Several formation modes can in principle reproduce the 
observed properties, but the old inferred stellar population 
ages restrict them to the first 4 Gyr. Unfortunately, the dif- 
ferences between plausible processes (e.g. major merging, 
satellite accretion or dissipational formation) blur at suffi- 
ciently high redshift - when gas absolutely dominates the 
baryonic content of (proto) galaxies - making the details of 
the formation process hard to identify. 

We can however put some constraints on the relevance 
of these main processes. For instance, we just pointed out 
that the major merger scenario fails to reproduce the ex- 
traordinarily rich GCS of NGC4594 - unless progenitors 
were fundamentally different from any surviving disc galaxy 
with measured GC specific frequency. Low-mass satellite ac- 
cretion cannot explain either the high metallicity of the 
outer spheroid or the extraordinarily populous metal-rich 
GCS. Both properties are however nicely consistent with a 
picture where cold gas flows dominate the first formation 
stages (Oser et al."'2010'). A two-stage dissipational collapse 
scenario (e.g., Forbes et al. 1997; Harris & Harris 2002) nat- 
urally results in an almost coeval, bimodal GCS with old 
ages and an MDF for the spheroid outer regions matching 
that of the metal-rich GC sub-population. Furthermore, re- 
cent hydrodynamical simulations suggest that this growth 
mode is preval ent for massive ga laxies inhabiting low-density 
environments (jOser et al.ll2010h . The formation of the em- 
bedded disc remains a puzzle, but whatever the detailed 
mechanism, it requires the infall of an enormous amount of 
high-angular momentum pre-enriched gas in order to repro- 
duce the sol ar and super-solar me tallicities of the outer and 
mner rmgs (jEmsellem et al.lll996l ). 

Finally, we stress that our study suggests that what 
has been called 'bulge' in NGC 4594 is likely not a 
classical bulge. We show indications that it is either a 
multi-component spheroid, or a system resembling regu- 
lar intermediate-mass elliptical galaxies. Ultra deep, multi- 
wavelength imaging of Sombrero could reveal new clues 
about its extended spheroidal component. 
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APPENDIX A: FITTING THE RINGS 

With the aim of testing whether the results and fits per- 
formed in Sect.[2]are affected by the fact that such fits do not 
include components to model the two main stell ar rings in 
Somb rero, we used the latest version of GALFIT (Peng et al. 
I2OIOI ). which allows one to model structures such as rings. 
The rings are described as elliptical Gaussian functions with 
inner and outer truncations. We first tried to start the fits 
from scratch. However, the parameters describing the rings, 
if left free, result in model components which are essentially 
merged with the galaxy disc, thus failing to reproduced the 
bumps observed in Sombrero's surface brightness radial pro- 
files, in particular along the disc major axis. We thus had 
to fix all parameters describing the rings, after many man- 
ual interactions. For this reason, and because the resulting 
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Figure Al. Left: surface brightness radial profiles (in arbitrary 
units) obtained through ellipse fits to the IRAC image and to the 
GALFIT model images of a decomposition including bulge, disc and 
two rings. Right: the same radial profiles, but obtained through a 
cut along the disc major axis. Lower panels show residual profiles 
after subtracting the full model profile from the galaxy profile. 
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Figure A2. Same as Fig. lAll but for a decomposition including 
bulge, disc, halo and two rings. 

models for bulge, disc and halo do not change significantly 
after including the rings, we keep using models bd and bdh 
from the analysis above as reference models. 

The results from these new fits are shown in Figs. lAll 
to IA3I The new model with bulge and disc plus rings is 
called BDR, while the corresponding model with an extra 
halo component is called bdhr. Figure lAll should be com- 
pared with Fig. [21 while Fig. IA2I should be compared with 
Fig.O As mentioned, the models for bulge, disc and halo do 
not change significantly after adding the ring components. 
The variations observed at each structural parameter are 
within the uncertainties found with budda for models bd 
and BDH, which are typically ^ 10%. 

Figure IA3I should be compared with Fig. [71 One sees 
that after the inclusion of the rings it is still evident that 
the extra halo component results in an improved fit. 
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Figure A3. Images of the full GALFIT models (top) and the corresponding residual images (bottom) after subtraction of the models 
from the IRAC image. Left: model bdr (bulge+disc+rings); right: model bdhr (bulge+disc+halo+rings). The images are in the same 
scale and display levels as in Fig. [7] 



